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Abstract 

We give the exact construction of Riemannian (or stringy) instantons, which are classical 
solutions of 2d Yang-Mills theories that interpolate between initial and final string configu- 
rations. They satisfy the Hitchin equations with special boundary conditions. For the case 
of U{2) gauge group those equations can be written as the sinh-Gordon equation with a 
delta function source. Using techniques of integrable theories based on the zero curvature 
conditions, we show that the solution is a condensate of an infinite number of one-solitons 
with the same topological charge and with all possible rapidities. 



1 Introduction 



In this paper we intend to prove the existence of Riemannian or stringy instantons. The name 
is due to the fact that these are classical solutions of a 2D U (N) YM theory that interpolate 
between initial and final string configurations. In other words they describe Riemann surfaces 
with punctures, where the latter represent the asymptotic entering and exiting strings. In the 
simplest case {N = 2), proving the existence of Riemannian instantons amounts to finding exact 
solutions of the sinh-Gordon equation with a delta-function source (see below for more details) . 
To be definite let a = a{z) be a polynomial in the complex variable z, with distinct roots, and 
let us introduce a new variable C defined by ^ = -y/a, where z = e^, w being the coordinate on 
an infinite cylinder. The equation we want to solve is 

d(;d^u - 2c/2sinh 2u = -^6{a){dt;a){d^a) . (1.1) 

which has to be understood in the sense of complex distribution theory. In this equation g is a 
constant coupling. An equivalent way to state (^]^) is to write the usual sinh-Gordon equation 

d(;d^u - 2c/2sinh 2u = (1.2) 

and to look for solutions which, near the zeroes of a, behave like 

n~— — ln|a|, for a~0. (1-3) 

We will also require that u vanishes as z ^ and z — > oo. This kind of equation was met for 
the first time in the framework of Matrix String Theory in |[l|, |2|, ^, ^, and solutions to these 
equations, satisfying ( |1.3| ), were shown to exist only numerically. In particular in [Q] this was 
done in the framework of a square lattice approximation, assuming a very simple form for a. 
The same type of solutions appeared in the context of form factors and correlation functions for 
the Ising model in [|[. 

This paper is devoted to proving the existence of solutions to the above equations, with the 
desired boundary conditions, in an analytic way, and to give their closed expressions in terms of 
the modified Bessel function Kq. In fact, the validity of the solution relies on some non- linear 
differential identites satified by integrals of Kq, which to our knowledge, have not appeared in 
the literature. 

The central ideas of the proof is (1) to use the Leznov-Saveliev approach and (2) to view the 
solution as a condensate of solitons. 



1 



More in detail, we begin by writing the sinh-Gordon equation (1.1) in terms of zero curvature 



conditions, which include besides the usual Lax-Zakharov-Shabat equation a second relation 
leading to non-local conservation laws. Such generalized zero curvature conditions follow from 
the ideas proposed in Q to study integrable theories in any dimension. Once we have the 
equations of motion written in terms of a zero curvature condition, we utilize the Leznov- 
Saveliev method to construct the corresponding Riemannian instanton solution. This method 
uses the fact that the dynamical variables of the system are contained in the zero curvature 
potentials: the sinh-Gordon ip field appears as a parameter of the group element we use in order 
to write the flat connetion as 

= -dfj,WW~^ = function of a group element 7, 7 = e'^^^'^"'^. (1.4) 

Thus, due to the path independence encoded in F^^, = 0, we are able to write the group element 
W in distinct forms. This, compounded with some properties of the Kac-Mooody algebra, leads 
us to a simple algebraic relation for the "group parameters" <p and u, i.e. 

(A I 7-1 I A) = (A I 7+Af+MrS- I A) (1.5) 

where the elements 7±, A'^^, and M_ have nice properties in terms of Kac-Moody algebra rep- 
resentations. To determine the two parameters ip and u, we make use of two highest weight 
representations. This provides us with a relation for ip in terms of the expected values of the 
group elements and M_. So, solving the sinh-Gordon equation is equivalent to furnishing 
these two elements. Once this is done, we conveniently choose the parameters in such a way 
that the boundary conditions (|1.3| ) are satisfied. 

The key point in this construction is that in order to obtain the desired solution we must 
choose the constant group elements of the solitonic specialization of the Leznov-Saveliev con- 
struction as an infinite product of exponentials of vertex operators. The product is in fact a 
continuous one, since it involves all possible values of the rapidities of the one-solitons. In ad- 
dition, all one-solitons entering the expansion have the same topological charge. This leads us 
to interpret such configuration as a condensate of solitons. As it was realized in this type of 
solution can be written as a Fredholm determinant, and our solution is similar to the one found 
in [||, where correlation functions of the Ising model were obtained in terms of r^^^ functions of 
the sinh-Gordon model, with N 00. In order to arrive at the true solution a continuum limit 
must be taken for the condensate of solitons and the Fredholm determinant must be rewritten as 
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an infinite series of integrals, whose convergence conditions are studied and particular normal- 
izations are fixed in order to satisfy the required boundary conditions. The solution heuristically 



derived in this way is finally shown to satisfy the equation (|l.l| ) or (1^) plus (|l.3| ). 

Before we enter into the very existence proof it is worth reviewing the framework where the 
eq. (|]|) arises and plays a fundamental role. Matrix String Theory (MST) |^ |9|, |l^ is the 



theory that arises upon compactifying Matrix Theory |11| on a circle, |12[| . It is expected to be 
a nonperturbative version of type II A string theory. An attempt to substantiate such conjecture 
was started in ||l|, ^ and completed in ^, where it was shown that a correspondence 
between MST and type IIA theory exists not only at the tree level, but that actually MST 
contains the full perturbative expansion of type IIA string theory. It was in this context that 



(|1.1| ) appeared. 

Looking for classical solutions of MST that preserve half supersymmetry, the following system 
of equations was found 

Dy,X = 0, Du,X = (1.6) 

where Fww denotes the curvature of a connection Au, and A^, while X is an x matrix and 
X its hermitean conjugate. Duj,Dw denotes the covariant derivatives with respect to Aw,Au,. 
(|1.6| ) may be called Hitchin equations, because they were discussed first by Hitchin in a different 
context [0], or Riemannian instanton equations because of their geometrical interpretation. To 
elucidate this terminology and the importance of these equations let us consider the simplest 
case, in which the gauge group is U{2). The problem to be solved is finding a couple {A,X) 
that satisfies ( |1.6| ). To this end we choose the following ansatz 

X = Y-'MY, A^ = id^Y\Y-^)^ , (1.7) 

where y is a suitable matrix G SL(2,C), and M is the following 2x2 matrix 

where a is a function on the complex plane. As a consequence of the equation D^X = 0, it 
follows that du,a = 0, i.e. a is holomorphic in z (at least for finite z). As explained above, we 
will assume that a is a polynomial in z with distinct roots. Now, given such an a we want to 
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find Y so that is satisfied. We parametrize y as y = ( ) wliere p = ^ + 4/n|a|, 

\ e~P 



and tt is a function to be determined. Then using (1.7) we find 

/ ae-2p\ /I 
^ = ( 2p ) ' = ( „ J (1-9) 




Now it is easy to verify that the first equation in (1^) imphes 



2d^d^p-gHe^P-\a\^e-^P) = (1.10) 



Inserting the explicit form of p and the change of variable w ^ s.t. 4^ = ^/a, one can rewrite 



(P^ as (O)- If 

M is a smooth solution of this equation, the couple {X, A) is a solution of (0) 



which is smooth everywhere except perhaps at infinity. Now, the important point is that the 
matrix M represents a branched covering of the z-plane. This is seen by diagonalizing M by 
means of a matrix in SL(2, C): 

M = SMS-\ M=[ , S=^[ , A. (1.11) 

The two eigenvalues of M represent the two branches of the equation 

y' = a (1.12) 

which is the defining equation of a hyperelliptic Riemann surface, with branch points corre- 
sponding to the roots of a. Therefore the solution at issue represents a Riemann surface, which 
justifies the adjective in the name Riemannian instanton. The instanton nature of this solution 
is due to the fact that ( |1.6| ) are the two dimensional reduction of the YM self-duality equation 
in 4D. In [jl5| this example was analyzed in great detail. It was shown there that, if the u 
solution of ( p..lj ) satisfies the boundary conditions stated at the beginning of this introduction, 
the matrix X, outside the branch points of a and when g ^ oo , becomes M up to a unitary 
transformation. This fact plays a crucial role in establishing the correspondence between MST 
and type IIA theory, see [^, 

It is therefore of upmost importance to establish the existence of the above solutions of ( [L.6| ) 



and therefore of the corresponding u solutions of (IT). On the other hand, it is clear that 
showing the existence of exact solutions of ( |1 . 1| ) is a problem interesting in itself. 

The paper is organized as follows. In section 2 we apply the Leznov-Saveliev method to 
our problem and define the soliton condensate. In section 3 we verify that what has been 
heuristically constructed in section 2 is in fact the looked for solution. A few appendices are 
devoted to clarify some technical points utilized in the course of the proof. 



2 Construction of solution through Leznov-SaveUev algebraic 
method 



In this section we review the Leznov-Savehev method [16] for construction of solutions of Affine- 
Toda type theories, based on the zero curvature formulation of two dimensional integrable 
systems. Even though we work in two dimensions we stick to the point of view of higher 
dimensional integrable models, which can be constructed from two potentials, A and B Q. 
Among other things such an approach leads to the construction of new conserved currents, not 
obtained from the usual two dimensional formalism 

2.1 Zero curvature condition 



The eq. (1.1) admits a representation in terms of the following zero curvature conditions 



F^,u = d^A,-d,A^ + [A^, A,]=0 (2.1) 
D^'B^ = d^'B^ + lA^', B^]=0 (2.2) 

In two dimensions the condition (|2.1[) is the well known Lax-Zakharov-Shabat equation. In 



dimensions higher than two the eqs. (|2.lD and (2^) were shown to be sufficient local conditions 



for the vanishing of the generalized zero curvature equations relevant for higher dimensional 
integrable theories Here we apply eqs. ( |2.1| )-( ^l2|) for the two dimensional model (1.1), and 



show that the relation (|2.2| ) leads to non-local conservation laws. The procedure applies equally 
well to a wide class of two dimensional integrable models, like the abelian and non-abelian 
Toda models (affine or not), possessing a representation in terms of Lax-Zakharov-Shabat eq. 



(2.1). That equation can be enriched by the extra conservation laws ( |2.2D without any further 
restriction in their dynamics, and we plan to analyse that in more details in a future publication. 

Let Q be an affine sl{2) Kac-Moody algebra, with basis H'^,T^,D,C, satisfying the com- 
mutation relations ( |B.113D . We take the local zero curvature potentials as 

= -9^77"^ + ^-1 A^ = 7^17"^ 

B^ = (E.i) ^ {lE,j-') (2.3) 

where 

7 = e^^°+^^ (2.4) 
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and 



El = gTl+ga{z) T}_ 



= ga{z) T-'+gT° 



(2.5) 



Following Ip we have taken Afj_ and to belong to a non-semisimple Lie algebra formed by the 
sl{2) Kac-Moody algebra and an abelian ideal which transforms under its adjoint representation 
P'^, i.e. 



rpm pip ('7^" 



pip (^m^ ^ pi, ^pn 







(2.6) 



With this choice, the zero curvature conditions (2.1) lead us to the system of differential 
equations 



g' (e'^- 



a |2 



g^ \a\^ 



(2.7) 
(2.8) 



The identification 



turns (|2.7D into the equation 



ip = 2p = u + — In \ a 



(2.9) 



(2.10) 



This is exactly the sinh-Gordon equation ( |1.1| ) once we make the change of variables w — > 
such that (see Appendix |^) 

^ -^ = ^/^ (2.11) 



dw 



As we point out in the introduction, the sinh-Gordon equation with source is equivalent to 
the homogeneous equation together with the following boundary conditions 



u ~ — In la 
2 ' 



Lp ~ finite 



for 



a ~ 



(2.12) 



So, u must diverge logarithmically at the zeroes of a, and from (2^) it follows that (p should be 
finite there. On the other hand, far away from any zero of a we need 



n ~ 



(f ^ — In \ a 



for 



a ~ 00 



(2.13) 
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One can check that the condition ( |2.2D is trivially satisfied by the potentials ( |2.3| ), i.e. it 
holds true for any field configuration including those which are not solutions of the equations of 
motion (|2.7| )- ( 12.81) . However, due to ( ^.l| ) it follows that the connection is flat, and so there 
is a group element W such that 

Af, = -d^WW~\ (2.14) 
Consequently, it follows from (^^) that the currents 

= W-^B^W (2.15) 



are conserved 



^^'J^ = (2.16) 



The group element W only exists for field configurations that satisfies the equations of motion, 
and it is non-local in the field variables. Consequently, the currents ( p. 15 ) are non-local. We 



intend to study the properties of these currents, in a wide class of models, in a future publication. 
2.2 Leznov-Saveliev construction 

The Leznov-Saveliev method uses some features of the Kac-Moody algebra representations, as 
well as some geometrical properties due to the zero curvature condition (|2.l|) . Since F^j^p = 0, it 



follows that the integration necessary for obtaining the group element W , introduced in ( 2.14| ), 
is path independent. That allows us to write this group element in different ways. We choose 
two different paths in order to integrate ( 2.14| ), and get 



W = gi=^g2 7 = 5192"^ (2-17) 

This means that we have two equivalent decomposition for the potentials A^: 

d^gig^^ = d^WW-^ = 0^77"' + 79m5252"S"' (2-18) 

Next we recall some aspects of the sl{2) Kac-Moody algebra. An affine Kac-Moody algebra 
possesses integer gradation such that |l^, |l8| 

Q = Q+ ® Qq ® Q- (2.19) 

where Qo^± are the eigensubspaces of the grading operator 

Q = ^H° + 2D; [Q,Gn]=nGn (2.20) 



with 

G+ = ^Gn G- =^Gn (2.21) 

n>0 n<0 

Observe that, in terms of the eigensubspaces Go,±, the potentials ( |2.3D have the decomposition 



Au, G ^4 



Bu, & P'^i G+ 



(2.22) 



We shall analyse how this potentials acts on a representation of the algebra. Consider a 
representation of the sl{2) Kac-Moody algebra. Inside it, there is a highest weight state | A) on 
which 

r_;'|A) = n>0 T!'|A) = n>0. (2.23) 

The scalar and spinor representations of the affine sl{2) Kac-Moody algebra have highest weight 
states I Ao ) and | Ai ) which satisfy 

C7|Ao) = |Ao) i7°|Ao) = (2.24) 

C I Ai ) = I Ai ) H° \ Xi) = \ Xi) (2.25) 

Therefore, using (|2T1) , (|2^ and ( ^) we see that the elements gi, §2 were taken in such 
a form that 

dingiQi^ I A ) = ( A I 8^0292^ = (2.26) 

This means that these group elements evaluated on a heighest weight representation are holo- 
morfic functions of the coordinates 

g^^\X)=i(w) {X\g2 = i{w). (2.27) 

With these ingredients at hand, we are interested in evaluating the expectation value 

(A I7-M A) = (A |52<7ri| A) (2.28) 

that descends from ( |2.17| ). To this end we perform the Gauss type decomposition 

gi = N-f-M- 52 = M7+iV+ (2.29) 

with 

7± G exp (Go) N, N+e exp (g+) M , £ exp (g^) ■ (2.30) 
So, using (|2l|), (|1|), (|239D and ([233| ) we compare the gradings of the potentials and of the 



group elements gi, g2, and with these we determine the action of 7 (2.28) on the representation 
states. 



Comparing the gradings of the potentials with the ones derived from (2.18) - ( |2.29| ), we 



get 



as weh as 






r-i 



disM_MZ^ = 
d^N+N-^ = 



d^NN~^ = 7^17-1 (2.31) 
dyjMNr^ = -7-^£^_i7 (2.32) 



So, as a consequence of ( |2.28| ) and ( p.29 ) the expectation values are 

(A 17-^1 A ) = ( A I 7+ (?D) iV+ {€)) MZ^ {w) 71^ {w)\ A ) 



(2.33) 



(2.34) 



This expectation values depends on two parameters 7±. A'^+ and M_ are determined from ( 2.33| ). 
Denoting 



7+ = exp {e+ [w) + e+ {w) C) 
one gets from (^) and (|2.34|) 



7- = exp {9- {w) + e- {w) C) 



(2.35) 



{M\N+ {€,) MZ^ {w) I Ao ) 



g-^-^ = e^+-e-+i+-i-^Xi\N+{w)MZ^{w)\ \i) 

At this point we are able to write the general solution of the model 

^ (Ai \N+{w)Mz\w)\ Xi) ^e^_e. 
(Ao I N+{w)Mz\w) I Ao) 
e"'^ = ( Ao I iV+ {w) MZ^ {w) \ Aq ) e«+-«- 



(2.36) 
(2.37) 



(2.38) 
(2.39) 



Notice that the group element W introduced in (2.14) have to be regular, since otherwise 
will not be flat. Remember that in order for A^ to satisfy the zero curvature condition it 
is necessary that the derivatives commute when acting on W . It then follows that the group 
elements gi, 52 and 7 have also to be regular. By regular we mean a quantity such that derivatives 
dw and dw commute when acting on it. From (see []T9| ) 

(2.40) 



5..-'=-^ = (-l)'=^5('='°) (.-,.) 



one observes that 



(2.41) 



Therefore, the fields and as weh as the parameters 6± and .^-i- can have log singularities, so 
that the corresponding group elements will have at most poles. 



The general solution of (^j), and so of is given by ( 2.38| ). Therefore the Riemannian 



instanton should correspond to some particular choice of the parameters of it. Due to (p.5|), 
(|2^ and (|2^ ) we have 



9^7V+iV+i = -5 (^e-2e+r0 +a(7l;)e2^+rlj (2.42) 
d^M.MZ^ = -g (a{w)e-^^-T^^ + e^^-T^) (2.43) 



This means that we have a solution to the model once we specify the parameters 6±, and 
solve (|1|), ( ra) for Af_, N+. 

2.3 The Riemannian instanton solutions 

As we have seen, the general solution of the model (|2.7|), or equivalently (pTiP, is given by ( p.38| ). 
We now have to choose the parameters and integration constants of the general solution, in order 
to obtain the Riemannian instantons solutions with the properties described in the introduction. 
We begin by choosing the functions 9± as 

6'+ = --lna 6i_ = -lna (2.44) 

4 4 

That choice simplifies the integration of the elements and M_. Indeed, ( |2.42| ) and ( |2.43| ) 
become 

d^N+N^^ = -gy^n{m) h d^M.MZ^ = -g./^ 6_i (2.45) 

The operators bi and 6_i are elements of a Heisenberg subalgebra of the sl{2) Kac-Moody 
algebra |17, 18 1. That is like an algebra of harmonic oscilators, i.e. they are generated by 

62n+l = Tl + ri+^ [62m+l , &2n+l] = C(2m + l)6m+n+l, (2.46) 

We can then integrate ( p.45| ) 



N+ = ^'h+ i+ = -g dw = -aC 



M_ = e^- ^-^h_ I- = -g / dw {w) = -gC (2.47) 



where we used the change of variables to ( |2.11 ), and where h± are constant group elements 



obtained by exponetiating the affine Kac-Moody algebra (integration constants). 
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We now return to the general solution ( |2.38| ). With our choice of 9± ( |2.44| ) and in view of 
dJ) we get 

(Ai I Af+M~i I Ai) , , 

e"" = ^ ~ ' 2.48 

(Ao I iV+MZ^ I Ao) 

Here we come to a crucial point in the construction of the Riemannian instanton solution. 
As it is well known |20, the one-soliton solutions are obtained by taking the integration 



constants h±, such that hj^h_^ = e^'^^\ where T^(^) is an element of the Kac-Moody algebra 
which is an eigenstate of the oscillators 62n+i) i-e. 



[62n+i, V{fi)] = -2^^"+^ (2.49) 

The construction of the operator V{^) is explained in appendix Its expression in terms of a 
special basis of the Kac-Moody algebra is given in ( 8.1291) . However, the nice properties of such 



operator are best appreciated in the principal vertex operator representation of the Kac-Moody 
algebra. Its form in that representation is given in ( B.151| ). An important relation satisfied by 



V{^) is given in ( B.158 ). From it we observe that 

V{^m)V{v) ^0 for (2.50) 

That implies that the exponential e^'^^^ truncates in first order, and so we do not have conver- 
gence problems in our expressions. Such property is what makes the vertex operator representa- 
tion to deserve the name of integral representation |]l^, |l8| . It also explain the truncation of the 
Hirota's expansion of the tau functions, since those are nothing more than special expectation 
values of V{fi) in the states of the vertex operator representation [pl| ]. 

If one takes a (z) to be constants and choose the integration constants h±, such that hj^hz} = 
e^'^^\ then one obtains from ( 2.4^ ) the one-soliton solution of the sinh-Gordon equation (by 



taking u ^ iu one gets the sine-Gordon one-soliton). The parameter ^ is related to the rapidity 
6 of the soliton through /i = ee^, with e = ±1. It is e what determines the sign of the topological 
charge (in the case of sine-Gordon) and what makes the difference between the soliton and 
anti-soliton solutions. 

The n-soliton solution is obtained by taking h^hZ^ as a product of those exponentials, i.e. 
h^hZ^ = niLi^^^'^*^- As we now explain the Riemannian instanton solution is obtained by 
taking h^hZ^ to be a continuous infinite product of exponentials e^^'^*^. In fact, we shall take 
the product in such a way that exponentials for smaller values of appear on the left, and we 
vary /ij continuosly from zero to -|-cxd. In addition, each value of /ij appears only once, without 
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repetition. So, what we have is an A^-sohton solution, with ^ oo, where aU the rapidities 
appear once, and we do not have a mixture of sohton and anti-sohtons since the /ij's are ah 
positive. Therefore, we have some sort of sohton condensat^. 

In order to build up the solution we start with an infinite discrete product of exponentials, 
and later take the continuous limit. So, we take the constant h+,h- in ( 2.47 ), ( ^.47 ) to be 



h+hZ^ = f[e 



(2.511 



From this point on we explore some useful properties of the vertex operators and its action on 
the representation states, which will allow us to end up with a closed expression for the solution 
u. So, once the constants hj^,h^ have been chosen as in (|2.51 ), the solution ( 2.48 ) depends on 



(A \N+MZ^\ A) = (A \ X{{l + e^^^^^V{tXi))\ X] 



i=l 



p{fi,) = 2g{fiU+ - ] (2.52) 



where we have used ( ^ ) and ( |230| ). Using ( pl60| ) and ( |Bl5^ ), we can expand ( ^ ) 

in terms of sums as 

2 



(Ao I N+MZ^ I Ao) 



1 + ^e/3(M.) + ^e^^'^^^e 



i<j 



IJ'k + IJ'j 



(Ai I N+MZ^ I Ai) = 1 - ^e^(^^) + ^e^^'^'^e 



i<j 



fXj + fit 



fJ-k + fJ-j 



+ ... 



(2.53) 



+ ... 



i<j<k 

Remember that we take the /Uj's to be real and positive, and that two //j's never coincide. 

Expressions (2.53) can be written in the form of Fredholm determinants. A similar result 
was found in ||5[, where exact correlation functions of the Ising Model were shown to be related 
to the tau functions of the sinh-Gordon model. Following H] we get 



(Ao I iV+Mr^ I Ao) 



(Ai I N+MZ^ I Ai) 



det (1 + W) 
det (1 - W) 



(2.54) 



^We are indebted to Olivier Babelon for pointing out to us that the Riemannian instanton solution should 
have such structure. His intuition came from his experience with D. Bernard on the calculations of form factors 
and correlation functions for the Ising model jsj 
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where W is the matrix 



Using (|2.4g| ) one then gets 



^ ^ e/3(M.)/2 \/?^g/3(M.)/2 (2.55) 



where we used In detM = Tr In M. Expanding the logarithm one gets 

Tr 



oo 



u = 2y^ (2.57) 

2.4 The continuous limit 

As we have said, we want to take the hmit where the infinite product of exponentials in ( |2.51| ) 
becomes a continuous one. In order to do that we take the label i of the parameter /ij to be the 
rapidity 9 of the soliton, and let it run from — oo to oo. We then have that ^ = . Then 



where A is a scaling factor of the integration measure. 
Then, from ( p. 55 ) it follows that^ 



(2.59) 



Therefore ( |2.57 ) becomes 

" = 2EyT^^2"+i (2-60) 



2n+ 1 

n=0 



with 



oo 



Jo /"I Jo m W + 1^2) [1^2 + IJ'3) if^N+f^l) 



In the case = 1 we have that 



2 Jo fJ- 



(2.62) 



where ATo is the modified Bessel function. In fact the above expression is valid for Re {gC) < 0. 



However, we shall take it to be valid for Re (gC) > too (see comments below ( p.68| )). 

■^Notice that the subindices of fi have a different meaning now. They label the variables giving the values of 
rows and columns of the matrices, and not the actual values of those as before. 
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Notice that the integrals 1^ are real. Indeed, from (2.52) we have 



(2.63) 



Therefore, one can undo the complex conjugation with the change of integration variables, 
l/^J-i, since ^ is left unchanged. In addition, /Xj/ (//j + fJ-j) — > fJ-j/ ifJ"i + fJ-j), and so 



the product of those terms in the integrand of ( 2.61 ) is left invariant. Consequently, the solution 
u given in (|2.60 ) is real (since A is real). 



We now want to analyse the boundary conditions satisfied by the solution (2.60). In order 
to do that we perform the change of integration variables 



In 



Mi 



l,2,...iV- 1 



(2.64) 



\i=l 



The integrals ( |2.61| ) become 
1 



•N-l 



cosh (i0i) cosh (i</)2) . . . cosh (|(/>Af-i) cosh ( ^ J2n=l 



where 



n 



(2.65) 



(2.66) 



1=1 



n=l 



n=l 



If Re {gC) < 0, the integral in u in ( |2.65| ) is the modified Bessel function Kq, and so one gets 



In 



2N-1 



')N-l- 



Ko{A\g\\C\ 



cosh cosh {^4>2) ■ ■ ■ cosh (i^jv-i) cosh ( i J2n=l 4>n 



where wn is given by 



WN = In {<P) In (- 



N-2N-1-1 j+l 
1=0 j=l i=j 



(2.67) 



(2.68) 



However, we shall take the expression ( 2.67] ) to be also valid for Re{gC) > 0. Such analytical 
continuation process will be justified later when we shall check the validity of the solution by 
directly replacing it into the equations of motion. Therefore, the solution (2.60) depends on 
and g through their norms only. 
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2.5 The boundary conditions 

For large arguments the modified Bessel function Kq have the following behaviour 

Ko (z) ~ e-' (^1 + 0^ ^ for large z (2.69) 

Consequently it is clear that 

In^O for ICHoo (2.70) 

and so, the u field does go to zero for large as required (see ( p.l3D ). 

The analysis for small C, is trickier. The reason is that taking | C | small does not guarantee 
that the argument of the Bessel function Kq is small, since wj\f can be infinitely large. However, 
in the region where | (" I w^TV diverges for small | C | the function Kq vanishes and so there is no 
contribution for the integral In- Therefore, we can use the following reasoning: Let | C | have a 
fixed infinitesimal value | (" |= e. We split the domain of integration in two regions, namely 

Do = region of (0i, . . . , (pN-i) where e4 | g' | ^/wn < 

Di = region of (0i, . . . , ^jv-i) where e4 | 51 | ^/wn > \/e (2-71) 
In the region Dq we use the fact that for small arguments, Kq diverges as 

Kq (z) ~ - In I (1 + O (z^)) for small z (2.72) 

and so 

r ^ f -In {2\g\\C\V^) 

= ttm-T / "01 • • • acpN-i ■ ■ ■ 



+ 



/ d(pi ■ ■ ■ d(pis[-i -7 : 

J Do cosh {\4>i) cosh {\4>2) ■ ■ ■ cosh [^(pN-i) cosh ( ^ X]n=l 

1 /• , , , , KQ{4\g\\C\^) 

/ d(pi . . . acpN^i -. 

cosh (^0i) cosh (^02) • • • cosh [^^'t'N-i) cosh ( \ Yln=l 



;2.73) 



Notice that in the region Di the argument of Kq never vanishes and so Kq is finite there. 
On the other hand, on Di we have 



and so — > 00 as e — > 0. But from ( 2.68 ) one observes that the only way for that to happens is 
that at least one the (/>j's should diverge. Therefore, the denominator of the integrand of ( p.73| ), 
in the Di region, diverges. So, one gets that the integral in Di in ( |2.73| ) vanishes for e ^ 0. 
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Consequently, the integral in Dq in (2.73) implies that 



In KAT In I CI for |C H (2.75) 

with 



1 dcpi... d(pN-i 



2^ ^ cosh (i0i) cosh (i(/)2) ••• cosh (i(/)Ar_i) cosh 



(2.76) 



Performing the integration using the fact that 



oo 



Kn = / dfpi . . . d(t)N iv = I dk 



°° Dill cosh 27r "^^ Hill cosh 



27r ./-oo '^^ Vi-oo ^"^coshf^ '''''' 



one gets 

tt'^"' (2n- 1)!! 



^2n^i = ^^^^ (2-78) 



Consequently we have from ( 2.6[1| ) that 



u -/(27rA)ln|C| for |CHO (2.79) 



vr 



where 



-^2"n!(2n + l) ""'"^^ 



(2.80) 



n=0 

Notice this series is convergent for < 1, and divergent for > 1. For = 1 the ratio test 
does not say anything, but one can check that it does converge there and / (1) = ir/2. Therefore, 
we must have | A |< l/27r. 

Our function a{z) is supposed to be a polynomial and to represent a hyperelliptic Riemann 
surface. As discussed in Appendix a good local coordinate near a branch point Zi is = 
^/z — Zi, i.e., near Zi we have z = Zi + . So, according to ( A.105| ), near the branching cut, we 
must have C, ~ a^/^. Therefore, 

u /(27rA)ln|o| (2.81) 

vr 



One can check that 



/|^)=f (2.82) 
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Consequently, in order to satisfy the boundary condition ( p. 12 ), we must set 



A = — (2.83) 

47r 



Therefore from ( p.60D , the desired solution to ( |l.l| ) is given by 

oo ^ 

u = 2y ^^"+^ , (2.84) 

^J2n+l)(2^)2"+i 

where l2n+i is given in ( 2.62| ) and ( p. 671 ). Consequently, it depends only on the combination 



4 I g II C I, and it is symmetric, on the plane (, under rotation around the origin C = 0. 

3 Check of the solution 

As we have seen in ( p.83| ), the value of the scaling factor A of the integration measure, introduced 
in ( p. 581 ), was fixed to l/Air. That was imposed by the behaviour of the solution at C = (or 
equivalently on the zeroes of a{z)). However, as we will see in this section, the solution holds 
true for any value of A, outside the zeroes of a{z). That is a consequence of very interesting 
non-linear differential equations satisfied by the integrals /2n+i defined in ( 2.67] ). Therefore, in 



order to emphasize those properties, we shall not fix A in this check of the solution. 

We begin by looking at the convergence of the series ( 2.60| ). As we have seen in ( |2.70| ), the 



integrals /2n+i go to zero for large arguments. In fact, from the numerical data in tables || and |2| 
of appendix one sees that they decay faster for larger index 2n + 1. Therefore, we should not 
have problems of convergence of the series ( p.60| ) for large arguments. For close to zero one 
can use ( p. 75 ) and (|2.78| ) and the ratio test to check that the series converges for A < ^/2/27r. 



Therefore, the series ( |2.84 ) for the final solution should converge everywhere. 

We begin by calculating the derivatives of the field u close to (" = 0. From ( |2.79| ) it follows 
that 

d^u ;^/(27rA) ^ (3.85) 



and so using (2.40) 



d^d^u^-f {271 A) 5 {C,C) (3.86) 

Some care must be taken here since we have been working with delta functions in different 
coordinate frames. In order avoid misunderstandings which can lead to inconsistencies in fixing 
A, we devote a discussion on this point in Appendix From ( [A.112| ) we see that 

d^d^u ~ -/ (27rA) 5 (C, C) = -/ (2^A) ^ 6 (a, a) d^a dfa , (3.87) 
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and using ( |2.82| ) and ( |2.83| ) we get 



TT 



(3.88) 



Therefore, we do have that ( |1.1D is satisfied at C = 0. 

Let us now evaluate the derivatives of u for (" 7^ 0. Since l2n+i depends on Q through the 
modified Bessel function Kq, we consider 

d-^d^K^ (4 U II CI V^) = (4 u D'ii^jv 5f I CI 5J CI K (4 U II C I ^/^) 

+ 4 I I ^d-^d(^ Id /iTo (4 U II C I V^) (3.89) 



Observe that 



and with help of ( 2.40 ) we obtain 

5c-5c I C 



' +7lClMC,C) 



4 I C I 4 

Since we are taking the point C = out, we can use the defining equation of Kq namely, 

z^K'^ {z) + zK'q {z) - z^Ko {z) = 

and (1^ and (|3^ ) to get 

d-^di^K^ (4 U II C I Vw^) = ^g^WN Ko (4 I 5 II C I Vw^) 



Therefore, from ( p.67D one has 

d-^d^iN = 45' Jn 

where 



Jn 



1 



3i . . . acpN-i- 



C/0 



w^Af -?^o (4 I 5 1 1 C I ^/w^) 



(3.90) 



(3.91) 



(3.92) 



(3.93) 



(3.94) 



cosh (i0i) cosh {\(p2) ■ ■ ■ cosh {^(pN-i) cosh ( i 0„ 



Notice that the relation ( 3.94| ) is also valid for = 1, with 

Ji=i^o(4UllCl) 



(3.95) 



(3.96) 



The reason is that from (2.62) we have that Ii = Kq (4 | (7 || C |), and so ( p. 93 ) becomes ( p. 94 
using the fact that 1111 = 1. 
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Using arguments similar to those leading to ( |3.93D , one can check that Ij^ and Jn satisfy 

z^I';^ (z) + zI'n {z) = z^Jn {z) (3.97) 

where z stands for the argument of those functions, i.e. 2; = 4 | 51 || |. 
We now have, from ( ^ ), (|3]94|) and (|3^ ), that 



d-^d^ u = 8g^Yl 2n + l -^'"^^ 'HC^Q f (2^^) (3-98) 

n=0 ^ 

where / (x) was defined in ( pO| ). 

So, replacing into equation ( |1.1[) , we get 

25M 4 ^ V+r " = -5 (C> C) / (27rA) - -5 (a) d^ad-^a (3.99) 

V n=0 / 

As we have already seen the vanishing of the r.h.s. of ( |3.99D fixes the value of A. Indeed, using 



( A. 112 ) we see we need to choose A such that / (27rA) = vr/G. But that is exactly what we had 



done in ( 2.82| ) and ( |2.83| ) to get the right boundary conditions. 



The l.h.s. of ( p. 99]) on the other hand vanishes for any A. That is an amazing result and 
involves special properties of the Bessel function Kq, or more precisely of In and Jn, which we 
could not find in the literature. Expanding the l.h.s. of ( |3.99| ) in powers of A we get 

A ^ Ko(4UII CI) =i^o(4UllC|) {h = Ji) 

A3 ^ J3 = /3 + 8/3 (3.100) 

- J5 = /5 + f/f/3 + f/f 

A^ - J7 = h + —lth + -hll + -llh + —ll 
' ' 9 ^■^ 9^ 5^ 45^ 



With help of ^M ) 

we get that the integrals /2n+i satisfy the following coupled non-linear 

differential equations 



^ X ^ 



T" -\- -T' 



h 
h 



X 



8lf 

40 2. ,32 5 
y /i /3 + y/i 



224 



56 



56 



1? Is H hl^ + —If 

9 1 ^ ^ 9 ^ 3 -r 5 1 



256 
45"^ 



(3.101) 
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So, the r.h.s. of these equations is what makes the difference between Kq (x) and the /2n+i's 
(see d^)). 

We did not find a way of checking relations ( |3.10l| ) analytically. In the appendix |C| we make 
a numerical check of them up to the equation for /y, and find that they are indeed true. 
Therefore, the configuration (^.84) is indeed a solution of (1.1). 
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A Hyperelliptic Riemann surfaces and delta function 

Let us consider the equation for an hyperelliptic Riemann surface 

y2 = a{z) = {z- zi){z - Z2)...{Z - z„) (A.102) 

There are branch points at z = zi, z = Zn- Let us suppose they are all distinct. Also the point 
at infinity is a branch point when n is odd. For simplicity let us suppose that n is even. There 
are two sheets, which are two copies of the z-plane. Now we draw a cut from zi to Z2, another 
from 2:3 to Zi and so on. The two sheets are now attached to each other through the cuts. By 
proceeding along a small circle around a branch point we will pass from one sheet to another 
and after Att we are back to the initial point. Passing to the string interpretation, it is evident 
that now we have a string of length 47r instead of 27r like initially. The string interpretation is 
as follows: we have initially two strings that interact successively n — 2 times and finally split 
into two separate strings as in the initial state. The Riemann surface we get in this way is an 
hyperelliptic one with two punctures representing the initial strings, two representing the final 
strings and (n — 2)/2 handles. Let us call it S. 

Let us return to (|A.102 ). y and z are coordinates of two complex planes, but, of course they 



can be considered as function over S. The coordinate z is not a good coordinate near a branch 
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point. A good local coordinate near a branch point Zi is ^, = -y/F^^Zi. I.e., near Zj we have 
z = Zi + z is not a good coordinate at infinity either, it must be replaced hy w = 1/ z. After 
these substitution we see that y is a meromorphic function, with n zeroes at the branch points 
and a pole of order n at z = oo on each sheet. 

Let us now consider the differential dz. dz ~ ^id^i near Zi, therefore dz has simple zeroes at 
the branch points. At infinity dz ~ w^'^dw, therefore it has a double pole there, on both sheets. 

Therefore the product ydz is a meromorphic one-form over S, with a single pole of order 
n + 2 at infinity. It makes sense to integrate this form along a path, and this is what we do 
when we write C = / y/adz, eq.(2.47). 

Let us now come to eq. ( |2.79D 

n~--/ln|C| (A.103) 

TT 



and to the behaviour 



u = -^ln|a| (A.104) 



which is required for consistency near a branch point. Near a generic branch point a ~ z — zq) 
a good coordinate is ^ = ^/a. In terms of this coordinate we have 

a ~ e^ C = j dzs/^ ~ j i'^di ~ e^ a ~ (^'^ (A.105) 

Therefore, in order that ( |A.103| ) be consistent with ( A.104| ) we must have / = f • 
Some attention must be payed to the definition of the delta functions, in order to avoid 
possible inconsistencies in fixing the value of /. Let us see this point in detail. 
Consider the good coordinate ^ and 

j d^e^c^lnlel = ly d\5{i,S,) (A.106) 

Since the contribution to the integral is only at the origin we can restrict it to the unit disk 
around the origin, and proceed in another way by applying Stokes theorem 

j d\ d^d^Mi\ = \j dH = ]^j>di-^ = \ j^J de = TTi (A.107) 

where the contour integral extends over the unit circle around the origin and ^ = e*^. 
If we repeat the same calculation with the 'bad' coordinate a we get 

j d^adada\n\a\ = ^ J d^a6ia,d) (A.108) 
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and 

J (fadadaln\a\ = ^ J (fada^ = ^ ^ da^ = ^ = 2m (A.109) 

The last steps are due to the fact that the angular integration for a extends over 47r since a ~ 
In a similar way for C, we will get 

j d''cdcd^\n\c\ = lj d''cd-^^ = ^fdc^ = ^l\e = 3m (A.llO) 

At this point it is judicious to make use of different symbols for these delta functions: S{^,^), 
which is the usual delta function, and Sa{a,d),S(^{C,,Q so that, roughly speaking, 

da{a,d) ~ 25{^,0, 5c(C,C) ~ m,0 (A.lll) 

In additon we must take into account the Jacobian factor due to the change of coordinates (a 
delta function transforms like the component of a one-form). In conclusion we have the relation 

SdCO = l6a{a,d)d^ad^d (A.112) 

B The afRne sl(2) Kac-Moody algebra 

The commutation relations of the affine sl{2) Kac-Moody algebra are given by 

= 2Cmdm+n,0 

[D,T^] = mT^ (B.113) 
In a highest weight representation we can take 

(if")^ = H-"" {Tl)^ = TZ^ (B.114) 

Let I A ) be a highest weight state. Then 

I A) =r:! I A) =il'' I A) =0 m>0; n>0 (B.115) 

Therefore, the norm of the state T^T™ | A), m > 0, is 

(A |r:!*T+"* I A) = (A I [T!",r+"^] I A) = (A l-ii-o + mCI A) (B.116) 
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Using 1 ^1^ ) and ( ^ ), one gets 



I Ai ) = 

Analogously, the norm of the state T^"^T^'^ \ A), m > 0, is 

^ ^ I ri.rnrj.mj.-mri.-m | ^ ) = 2( A | + mC) {-H^ + mC) - {-H^ + mC) I A ; 

Therefore, from ( \2.24) 



(B.117) 



T7'T7' I Ao ) = 



Using similar arguments one gets that 



I An ) = 



T°T° I Ai ) = 



(B.118) 



(B.119) 



For the applications we make in this paper it is useful to work with following basis for the affine 
sl{2) Kac-Moody algebra 



^2m+l 


= T"^ 


_j_ rpm+1 


(B.120) 


F2m+1 


rpm 

— + 


rj.m+1 


(B.121) 


F2m 




— -C^m,0 


(B.122) 



in addition to the generators C (central term) and D . The commutation relations are 



[^2m+l , ^2n+l] 
[^2m+l , F2n+l] 
[^2m+l , F2n\ 
[F2m+1 , F2n+l] 
[F2m+1 , F2n\ 
[F2m , F2n] 



C{2m + l)5m+n+lfl 

— 2F2(m+n+l) 

— '^F2(m+n)+l 
-C{2m + l)5m+n+lfl 

— 262(m+n)+l 



C2m5m+n,0 

The elements which diagonalize the adjoint action of the oscillators 62m+i ar^ 

oo 
n=— oo 

Indeed 

[h2m+l,v{^l)] = -2^?^+^v{^l) 



(B.123) 
(B.124) 
(B.125) 
(B.126) 
(B.127) 
(B.128) 



(B.129) 



(B.130) 



^We choose the normahzation factor to be —2 because we want the expectation values in (B.131) to be ±1. 
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Notice that the expectation value of such operator on a highest weight state gets contribution 
only from its zero mode. Indeed, using and (I^ID, one gets 



(Ao Ao) = l (Ai Ai) = -1 (B.131) 

B.l Vertex Operator Construction 

We now discuss the principal vertex operator representation of the sl{2) Kac-Moody algebra. 
It differs from the homogeneous vertex representation, i.e. the Fubini-Veneziano vertex, in the 
sense that the oscillators have only odd indices and so do not have zero modes. We introduce 

°° / -N N \ 

^(^) = E(i^^n-^6_n) N = 2n + 1 (B.132) 

n=0 ^ ^ 

oo 

P(z) = -Y^{z'^b^ + z^b^^) (B.133) 

n=0 

where the oscillators satisfy the commutation relations ( B.123| ). We then have 



w z 

P{z)Q{w) = :: P{z)Q{w) :: + ^ C \z\>\w\ 

Q{w)P{z) = :: Q{w)P{z) :: - ^ C \w\>\z\ (B.134) 

— z^ 

We now introduce the vertex operator 

V(z) = :: e"^(") :: (B.135) 

and so 

OL wz 

P{z)V{w) = :: P{z)V{w) :: + ^ CV{w) | ^ |>| u; | 

— 

rv unz 

V{w)P{z) =::V{w)P{z) :: . ^ CV{w) \ w \>\ z \ (B.136) 

We expand V{z) in modes as 

/rl7 7" 
^V(z) (B.137) 

That means that 6M^n and ^n^M can be written as the doble integral of the same term 





fdz z^ 


27riC J 


2711 z 


dCC 


fdz z^ 


2m C J 


2m z 



\A>\Q 



a (z 



AnbM = -{<Pi±h<P^TZ—]^^^^^[--P(^)^iO--+^2^C^(0} (B.138) 
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where we use 



Notice that 



dz 
27ri z 



P{z) = -6m M odd 



dCC" fdzz 



27ri C, J 27ri z 



\A>K\ 



dC C f dz z^ ' 
2Tri( J 2TTiz y |^|^|^| 



(B.139) 
(B.140) 











© 








© 





The only contribution comes from the poles on z = zt^; we eliminate the poles on z = 0. So we 
have 

dec f dz z^ ( aCz 



2.^C J 2mz ^::^mC)::+^CV(C) 



fB.Ml) 



The first term in r.h.s. is normal ordered, what guarantees that there is no contribution of this 
part on z = C- The second term has a simple pole on z = ib^. So, by the residue theorem we 
have 

- -|^/lf^C-V(c).|c/|^,-crv,c) 

= -aCAn+M (B.142) 



once M is odd and we used (lBl37|) , (^l39D . This relation corresponds to ( pl2l , ( pl25|) , but 
in order to reproduce the algebra, we must impose 



aC = 2 



[6m, An] = —2An+-M 



(B.143) 



Now we must verify the remaining commutation relations, namely [^mj^n]- To this, we must 
evaluate 



V(z)V(C) = :: V(z)V(C) :: e"^^ '""^S?) \ z \>\ C 



(B.144) 



from where we conclude that 

[Ami ^n] = 



dz / 


dC C 


27ri z J 


2iriC 






(z-C 


y 







:: V(z)V(C) :: x 



l^l>ICI 



■ c 



Kl>kl. 



(B.145) 
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Since we need double poles, we set ^ C = 2. But from ( B.143| ) we have aC = 2, and so we get 
that a = 2 and C = 1. Then 



dz 

27ri z 



2TTiC 



:: V{z)ViC) :: 



z-C 
z + C 



(B.146) 



where we take the contour to be 
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According with the Cauchy's integral formula 

1 f f{z) 



fizo) 



2m 



{z - zq)'- 



dz 



what implies 



dC 

Zm Q 



(B.147) 



(B.148) 



(B.149) 



where we have used the fact that :: V(— C)V(C) - = 1- Using ( B.139| ) one gets 

[^2m+l, -42„+l] = -4 (2m + 1) (5m+n+l,0 
[^2m+l, ^2n] = " 8^2(m+n+)+l 

Therefore, from (|B.143D and (|B.149|) , we see that zizAn/2 satisfy the same algebra as Fn, 
given in (B.124)-( B.128D . So we have reproduced the commutation relations for the sl{2) Kac- 
Moody algebra for C = 1. The arbitrariness in the sign of An is related to the normalization of 
the expectation values of the operators on the highest weight states. Since the operator V (z) 
introduced in ( p.l35 ) is normal ordered, it satisfy 

(A \V{z)\ A) = 1 (B.150) 



Then, in order to reproduce ( B.131 ) we must normalize V (z) differently in different highest 
weight representations. In fact, we have to make the identification 

V{z) = cxV{z) (B.151) 

where 

CAo = 1 CA, = -1 (B.152) 
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B.2 Vertex Operator properties 

Using (^l32D 



n=0 



n=0 



N = 2n + 1 



v(z)VH 



g2Q<(^) g2Q>(^) g2Q<(t«) g-2Q>{2) ^2 Q>{^) ^2 Q> 



g2Q>(^) g2Q<(u>)g-2Q>(^) 



g2Q>(^) 2Q<(to) e-2«>(^) 



g2Q<{«,)+[2Q>{^),2Q<(«,)] 
,2Q<M -2 ln(f±f) 



where we used 



[2Q>(z), 2Q<H] = -4 ^ 



°° — N M 

Z W 



n,m>0 



N M 



-M 



°° -N M 

Z W 



N M ' 

n,m>0 



n>0 



-2 In 



z + w 



z — w 



Then, from ( plsHD we have 



,2Q>iz) 2Q<{w) ^ 2Q<{w) p2 g> (^) -2 ln(f+f ) 



So, returning to ( [B.154 ) and using ( |B.157 ) 

V{z)Viw) = e2Q<(^) e2Q>(^) g2Q>H g-2 in(f±f ) 

2 



= :: V{z)V{w) :: 

By the same argument, 

V{z)V{w)V{k) = :: V(z)VHV(A;) :: 



z — w 



z + w 



\2/ ;\2/ ,\2 
Z — W\ I Z — k\ / W — k 



z + wl \z + kj \w + k 
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and so on. Therefore, from ( B.151| ) we get that 

(A I Vin,)...Vif,n) I \)=cl n ^^V^ (B.160) 



j>t=i ^ 

where we have used that fact that 

(A I :: V(/ii)...V(Mn) I A) = 1 (B.161) 

C The numerical check of equations ( |3.1U1| ) 

As we have seen in section ^, the existence of the exact solution rehes on the vaUdity of the 
non-Hnear differential equations (3.101) satisfied by the quantities /2n+i- We could not find in 



the literature any reference to those equations, and we were not able to prove them analitically. 
Therefore, we present here a numerical check of the first three equations in ( [j.lOll ), at some 
specific values of the argument of those functions. 

In order to perform the numerical calculation it is better to make a change of integration 
variables in the expressions for and J^. We define 

0. = lnii^ i = l,2,...N -1 (C.162) 



I - Xi 



Then 



We have that 



In addition we have 



Xi = —1 — > 4>i = —oo 
Xi = 1 4>i = CO (C.163) 



^ dxi (C.164) 



l-xi 




1-xf 



J2 , X j2 



2 mn (1 - 
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(C.165) 



We then get that 1^ and Jj^ given in ( p.67 ) and ( 3.95| ) respectively, become 



where 



VN 

1 



J-1 VN 

/=o j=i \i=j ^ i=j y ' ^' • 



and 



^ /N-l N-i \ 

= 2 n (1 + + n (1 - (^-168) 

\ j=l i=l / 

In order to evaluate the multidimensional integrals we used the Monte Carlo method. The 
program was written in Mathematica and did not make use of the command NIntegrate, 
since we found some limitations on it regarding the number of iterations. We found it bet- 
ter to write a routine to do it directly, using Mathematica random number generator and its 
package for Bessel functions. The functions In and Jn are defined as {N — l)-dimensional in- 
tegrals. We have evaluated them for = 3,5, and?, by sampling the integration region with 
10^ points, irrespective of the dimensionality of the integrals. We evaluated those functions for 
the arguments z = 0.001,0.01,0.1,0.5,1,5,10, and 100. In order to estimate the errors in the 
calculations, we evaluated each integral for each value of the argument, 11 times with the 10^ 
point sampling. We then calculated the average {v = ^'i=iVi/n) and the standard deviation 
(^' = Er=i {v^-vf/{n-l)). 

Using (|3.97| ) one can write the equations (|3.101 ) as 



where 



J2n+1 — L2n+l (C.169) 



Li = h + Ml 



Therefore, we check the relation ( C.169| ) numerically. The errors in L2n+i are given by 
ALs = A/3 

40 n 

ALs = Ah + —I^Ah 

AL, ^ Aly + ^lfAI^ + ^hhAh + ^lfAh (C.171) 
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z 


h = Kq 


0.001 


7.02369 


0.01 


4.72124 


0.1 


2.42707 


0.5 


0.924419 


1 


0.421024 


5 


3.6911 10-3 


10 


1.77801 10"^ 


100 


4.65663 10-"^^ 



Table 1: The numerical values of Ii for some values of the argument 



Since Ii = Kq, and since it is evaluated with the Mathematica package for Bessel function, its 
error was taken to be negligible as compared to the other ones. Its values, for the arguments 
used in the simulations are given in table ||. In table ^ we show the results of the numerical 
calculations of the integrals I3, and I/, with the corresponding standard deviations. 

The check of the relations ( C.169| ), or equivalently (3.101), is given in tables g I and H 
Notice that the agreement is quite good. Except for the values of J5 and L5 at z = 0.001, which 
differ by 1.45 standard deviations, all the other values agree well inside one standard deviation. 
So, our numerical calculations strongly indicate that the relations ( p.l01| ) should hold true. It 
would be very interesting to find an analytical proof of those equations. 
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z 


o 




i 


0.001 


26.237 ± 0.026 


155.52 it 0.44 


1027.0 ± 3.2 


0.01 


14.92 ± 0.01 


73.30 ±0.12 


390.4 ± 1.0 


0.1 


4.455 ± 0.002 


11.429 ±0.011 


30.118 ± 0.035 


0.5 


0.43063 ± 0.00044 


0.24327 ± 0.00024 


0.13817 ±0.00031 


1 


0.049517 ± 0.000035 


0.006656 ± 0.000010 


(8.974 ± 0.028) 10"^ 


5 


(4.4348 ± 0.0068) 10"^ 


(5.548 ± 0.018) 10-13 


(6.9 ±0.1) 10-1^ 


10 


(5.25 ±0.02) 10-15 


(1.589 ±0.017) 10-24 


(4.79 ±0.21) 10-34 


100 


(9.990 ± 0.065) lO-i^'' 


(2.21 ± 0.27) 10-222 


(3.8 ±2.8) 10-311 



Table 2: The numerical values of I3, I5 and I7 with corresponding standard 
deviations 



z 


h 




0.001 


2783. ± 338. 


2798.19 ± 0.03 


0.01 


863.9 ± 27.8 


856.82 ± 0.01 


0.1 


118.83 ±0.27 


118.831 ± 0.002 


0.5 


6.7521 ± 0.0074 


6.75033 ± 0.00044 


1 


0.64645 ± 0.00040 


0.646569 ± 0.000035 


5 


(4.4631 ± 0.0063) lO""^ 


(4.4665410 ± 0.00068) IQ-'^ 


10 


(5.022 ±0.021) 10^1'' 


(5.0220 ± 0.0022) lO-i^ 


100 


(9.051 ± 0.059) 10-133 


(9.0770 ± 0.0065) 10-133 



Table 3: Comparison of the numerical values of J3, and L3 with correspon- 
ding standard deviations 
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z 


^5 


^5 


0.001 


181810. ± 12383. 


199742. ± 18. 


0.01 


29616. ± 318. 


29529.2 ± 3.1 


0.1 


1260.3 ± 1.7 


1259.64 ± 0.17 


0.5 


12.35 ±0.01 


12.3505 ± 0.0052 


1 


0.26474 ± 0.00032 


0.264801 ± 0.000093 


5 


(1.5883 ± 0.0052) 10"^^ 


(1.5919 ± 0.0014) 10-^^ 


10 


(4.273 ± 0.046) lO'^^ 


(4.268 ±0.011) 10-23 


100 


(5.56 ± 0.67) 10-221 


(5.45 ± 0.46) 10-221 



Table 4: Comparison of the numerical values of J5, and L5 with correspon- 
ding standard deviations 



z 




L7 


0.001 


(6.40 ± .65) 10^ 


(6.5034 ± 0.0019) 10^ 


0.01 


505613. ± 4987. 


507197. + 166. 


0.1 


7753. ± 11. 


7753.3 ± 2.9 


0.5 


14.641 ± 0.027 


14.642 ± 0.013 


1 


0.07265 ± 0.00018 


0.072600 ± 0.000059 


5 


(3.932 ± 0.053) lO^^^ 


(3.9474 ± 0.0083) lO^^^ 


10 


(2.54 ±0.11) 10-32 


(2.542 ± 0.016) 10^32 


100 


(1.9 ± 1.4) 10-309 


(2.3 ± 1.0) 10-309 



Table 5: Comparison of the numerical values of J7, and L7 with correspon- 
ding standard deviations 



32 



References 



[1] T. Wynter, Gauge fields and interactions in matrix string theory Phys.Lett. B415 (1997) 
349 [hep-th/9709029]. 

[2] S.B. Giddings, F. Hacquebord, H. Verlinde, High Energy Scattering of D-pair Creation in 
Matrix String T/ieori/ Nucl.Phys. B537 (1999) 260 [hep-th/9804121]. 

[3] G. Bonelli, L. Bonora and F. Nesti, Matrix string theory, 2D instantons and affine Toda 
field theory, Phys.Lett. B435 (1998) 303 [hep-th/9805071]. 

[4] G. Bonelli, L. Bonora and F. Nesti, String Interactions from Matrix String Theory, 
Nucl.Phys. B538 (1999) 100 [hep-th/9807232]. 

[5] O.Babelon, D.Bernard, Int.J.Mod.Phys. A 8 (1993) 507 [hep-th/9206003]. 



[6] O. Alvarez, L.A. Ferreira, J. Sanchez Guillen, Nucl. Phys. B529(1998) 689, hep-th/9710T47 



[7] A.C.Scott, F.Y.F.Chu, D.W. McLaughlin, Proc. IEEE 61 (1973) 1443, L.D. Faddeev, L.A. 
Takhtajan, Hamiltonian Methods in the theory of Sohtons (Springer, Berlin, 1986). 

[8] L. Motl, Proposals on Nonperturhative Superstring Interactions, [hep-th/9701025]. 

[9] T. Banks and N. Seiberg, Strings from Matrices, Nucl.Phys. B 497 (1997) 41 [hep- 
th/9702187]. 

[10] R. Dijkgraaf, E. Verlinde, H. Verlinde, Matrix String Theory, Nucl.Phys. B500 (1997) 43 
[hep-th/9703030]. 

[11] T. Banks, W. Fischler, S.H. Shenker and L. Susskind, M Theory As A Matrix Model: A 
Conjecture, Phys.Rev.D 55 (1997) 5112 [hep-th/9610043]. 

[12] W. Taylor, D-brane Field Theory on Compact Spaces, Phys.Lett. B394 (1997) 283 [hep- 
th/9611042]. 

[13] G. Bonelli, L. Bonora, F. Nesti and A.Tomasiello, Matrix String Theory and its Moduli 
space Nucl.Phys. B554 (1999) 103 [hep-th/9901093]. 

[14] N.J. Hitchin, The self-duality equations on a Riemann surface, Proc. London Math. Soc. 
55 (1987) 59; Lie groups and Teichmuller space, Topology 31 (1992) 449. 

33 



[15] L. Bonora, Yang-Mills Theory and Matrix String Theory, Quantum Field Theory. A 20th 
Century Profile. A.N.Mitra ed. Delhi, 2000. 

[16] A.N. Leznov, M.V. Saveliev,Lett.Math.Phys. 3 (1979) 489; Group- Theoretical Methods for 
Integration of Non-Linear Dynamical Systems, (Birkhauser, Basel, 1992) 

[17] V.G. Kac Infinite Dimensional Lie Algebras, Cambridge UK: Univ. Pr. (1990). 

[18] D. Olive, P. Goddard, Kac-Moody nad Virasoro algebras in relation to quantum physics 
Int.Jour.Mod.Phys. Al (1986). 

[19] I.M. Guelfand, M.I.Graev, N.Ya. Vilenkin, Les Distributions, Tome 5, page 499. 

[20] D. I. Olive, J. W. R. Underwood and M. V. Saveliev, Phys. Lett. 311B (1993) 117-122, 
|hep-th/9212123 . 



[21] L. A. Ferreira, J. Luis Miramontes and J. Sanchez Guillen, J. Math. Phys. 38 (1997) 



882-901, |hep-th/9606066 



34 



